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Abstract 

The Leibler theory [L. Leibler, Macromolecules 13 1602 (1980)] for mi- 
crophase separation in AB block copolymer melts is generalized for systems 
with arbitrary topology of molecules. A diagrammatic technique for calcula- 
tion of the monomeric correlation functions is developed. The free energies 
of various mesophases are calculated within the second-harmonic approxima- 
tion. Model highly-branched tree-like structures are considered as an example 
and their phase diagrams are obtained. The topology of molecules is found 
to influence the spinodal temperature and asymmetry of the phase diagrams, 
but not the types of phases and their order. We suggest that all model AB 
block-copolymer systems will exhibit the typical phase behaviour. 



I. INTRODUCTION 



One of the biggest challenges for theoretical polymer science is to predict the phase 
behaviour of a polymer system with given molecular structure. The required theory should 
explicitely contain molecular parameters and be able to predict thermodynamically stable 
phases (whatever they are) and temperatures of transitions among them. Apart from the 
theoretical interest this problem is of great industrial importance. Since the mechanical 
properties of the phases are studied quite well it is useful to predict the conditions at which 
these phases are stable. In the absence of such predictions the experimental search could 
take even years. In the present situation of a rapidly increasing number of novel materials 
even a crude theoretical model could help a great deal in direct applications. 

Among the systems of interest are solutions of homopolymers, block copolymer melts, 
blends and so on. In this article we focus our attention on the A-B block copolymer melts. 
Such systems can only differ in the total content of one component, the length of molecules 
and the way the A and B blocks are arranged in one molecule. As a function of these pa- 
rameters and temperature, the following phases were observed experimentally: the lamellar 
phase (LAM), the hexagonal phase (HEX), the body-centered cubic phase (BCC) and the 
bicontinuous gyroid phase (GYR) However, the structure of these phases strongly de- 

pends on temperature. Just below the order-disorder transition temperature the chains are 
slightly stretched and the mesophases consist of periodic patterns with only a small increase 
in the content of one component. This situation is called the weak-segregation regime. On 
contrary, the strong-segregation regime (very low temperatures) is characterized by domains 
which are almost pure in one component. Both regimes are limiting cases and therefore can 
be studied analytically. In the intermediate regime there is no small parameter and only 
numerical investigations are possible 

The first weak-segregation theory of AB melts is by Leibler Q who considered a model 
Gaussian diblock copolymer. He was able to estimate the order-disorder transition tem- 
perature and predict the following sequence of transitions: DIS^BCC — >HEX— >LAM with 
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decrease of temperature. Later, after the GYR phase was observed experimentally 
Milner and Olmsted 0] used the second-harmonic approximation and predicted the GYR 
phase to be stable in between the HEX and LAM phases. The same conclusion was drawn 
in and in the Landau-Brazovskii theory of Hamley and Podneks ■ 

The method developed by Leibler was applied to other topologies of block copolymer 
molecules. Thus, Mayes and de la Cruz considered an ABA triblock molecule |ll|], de la 
Cruz and Sanchez studied star- and simple graft-copolymers ||12|; Dobrynin and Erukhi- 
movich developed a computer-aided tool to calculate phase diagrams for various topologies 
||13|| . However, in all these studies the authors started with a fixed molecular topology. To 
consider another topology one needed to go through the whole calculation procedure again, 
from the very beginning. In this work we show how the Leibler approach can be used for an 
arbitrary topology of molecules. Section I repeats some main points of Leibler's paper ^ and 
introduces a diagrammatic technique that helps to deal with molecules of arbitrary topol- 



ogy. The first step in this direction was made by Read |[T^, who invented a way to calculate 
structure factors of polymers of arbitrary complexity. We go further and show how to calcu- 
late higher vertex functions for such systems. Section II describes the way to calculate the 
free energies of different mesophases and discusses some of the approximations. We consider 
the face-centered cubic phase (FCC), the ordered bicontinuous double- diamond phase, the 
simple cubic (SC) and the square lattice (SL) phases in addition to the classical phases. In 
the Results and Discussion we apply the developed technique to a model dendrimer system 
which has recently attracted wide experimental [|15[ and theoretical attention. The phase 
diagrams of various dendrimers are presented and their features are discussed. The results 
for many more systems that were not included in this publication are presented as sup- 
plementary materials on the Internet: fittp : / /rug'md4 ■ chem. rug. nl/ ~morozov /research. htm^ 
Finally, we summarize our work in the Conclusion and propose some possible developments 
of the method. 
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II. MEAN-FIELD THEORY 



A mean-field theory of incompressible AB block copolymer melts was first formulated 
by Leibler [Q. In his approach a local deviation of the density of A component from its 
average value plays the role of an order parameter ilj{r). He found that the coefficients in 
the phenomenological free energy expansion 
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(2.1) 



can be connected with microscopic properties of the system. This possibility is based on the 
assumption that in a melt a single chain behaves as almost ideal. This constitutes the so- 
called random phase approximation (RPA) which expresses the monomer density correlation 
functions in terms of the correlation functions of a Gaussian chain. For the vertex functions 



in ( |2.1|) it gives 



r2(qi,q2) = 5(qi + q2) 



r3(qi, q2, q3) = -Gfh(\i, q2, q3) Sj{cii) - Sj(qi 



SaaM + 'S'BB(qi) + 2Sab{<^i, 
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2X 



X 5-i(q2) - 57i(q2)] [^,-i(q3) - S^'^^) _ 
r4(qi, q2, q3, q4) = lijki{<^i, q2, q3, q4) [SIaM - [Sjl{q,2) - SjB{ci2] 
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where 



lijkMi, q2, q3, q4) = -G'Sjtz(qi, q2, q3, q4) 



+ / S. 



-1 



G 



(3) 



(qi, q2, p)G'S(-p, q3, q4) + G'Sfc!„(qi, qs, p)G'iJi(-p, q2, q4) 

+ <-^i?m(qi5 q45 p)G'ifcj(-P) q3, q2) 



(2.3) 



The summation over repeating indices is assumed. The second, the third and the fourth order 
correlation functions of the Gaussian chain are denoted by Sij, OfX and CfXi, respectively 
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{i,j,k,l = A,B); S^j^ are the elements of the inverse matrix \ \Sij\\ ^. The Flory-Huggins 
parameter describing the strength of interaction between A and B monomers is denoted by 

X- 

A general method of calculating the Gaussian correlation functions was developed in 
Appendix B of |^ . It is based on the fact that the Gaussian distribution is fully determined 
by its first two moments. As a result the probability to find several monomers at given 
positions on the chain is expressed as a product of the 2-point correlation functions of the 
consequently neighbouring monomers along the chain. For example, 

Ti r2 r3 
> . • = P(ri, r2)P(r2, rg) 



The correlation functions Sij, G^jl and G^jli are calculated by averaging the Fourier 
transforms of expressions like one in the previous example over all possible positions of 
monomers along the chain. At this point one needs to introduce some topology of the 
polymer molecule in order to perform such an averaging. In |Q this was done for a diblock 
molecule. A triblock molecule was considered in ||lT|. More sophisticated architectures were 



studied in [|13|. We generalize this procedure for arbitrary topology of a polymer molecule. 

In the following we use the graphical method formulated in Appendix A. A correlation 
function of a given order is calculated by drawing all possible (but topologically distinct) 
diagrams of this order. Thus, 



This is just a schematic representation which shows that an n-th order diagram can con- 
tain up to n blocks. The corresponding contributions to the correlation function are obtained 
by connecting these blocks in all possible ways. The topologically different connections are 
presented in Appendix A. The contributions to the correlation functions are calculated in 
Appendix B. 

We want to note that not all diagrams contribute to the correlation functions of the 
system with given topology. One restriction is connected with the type of the correlation 
function we want to calculate. For example, the equation for Sab only contains diagrams 
with the first label placed on some A block and the second label placed on a B block. 
Therefore the first diagram in the general equation for Sij does not contribute to Sab- 
The other factor reducing the number of diagrams is a molecular architecture. Thus, for a 
linear multiblock molecules all 'star' and 'fork' diagrams are impossible and should not be 
considered. 

The final expressions for the correlation functions are obtained in the following way. 
Suppose we want to calculate G'ii,i2,...- To do this we need to place the first label on any 
block of the type ii, the second label on any block of the type 12 and so on. The resulting 
diagram will be one of the diagrams calculated in Appendix B. The corresponding expres- 
sion is written down. Then we remove all labels and place them again avoiding previous 
configurations. This procedure is repeated until we covered all possible arrangements of the 
labels. The obtained expression is the equation for the required Gji^ia,.... 

Once we have calculated the correlation functions for a given system, its free energy can 
be obtained with the help of ( p.l[ - p.3|) . In the next section this Landau-Ginzburg free energy 
will be a starting point for our analysis of the phase transitions. 

It is obvious that the formulated approach is general and within this framework all 



systems can be considered. 



III. CALCULATION OF PHASE DIAGRAM 

In this section we show how to calculate a phase diagram of a system once its free energy 
functional ( |2.1| ) is known. Apart from some technical details this was developed elsewhere 
™jT8|J7|OJTl-01P,|T0i. 



In the beginning we need to determine the order-disorder transition (ODT) temperature. 
The disordered phase can be described by the Landau-Ginzburg Hamiltonian with only a 
quadratic term left The probability to create a spatial density fluctuation ipi^) is 

given by 

1 



exp 



"2 / ^2(?)^(q)^(-q) 



(3.1) 



The typical life-time of such a fluctuation is r 



1 



r2(g) 



. At some g*, temperature (or 

Flory-Huggins parameter) and total content of A monomers in the molecule ^/ = ^ 
1^2(9*) reaches zero. The probability to create such a mode tends to unity, while its life-time 
Tg, diverges. It means that at this temperature the disordered phase becomes absolutely 
unstable with respect to the appearance of an ordered mesophase with typical domain size 
~ q^^. Therefore, the set of equations: 



r2(g*,(xA^). 







(3.2) 



dT2{q, {xN)s) 



dq 







d^T^iq, (xiV). 



dq"^ 



> 



determines the spinodal temperature (which is very close to the ODT value [^) and the 
typical lengthscale q^,. 

As the next step we choose a set of mesophases which will possibly enter the phase 
diagram. Then, for each mesophase we represent the order parameter fleld ip{c[) in such 
a way that it possesses the symmetry of the mesophase. This representation necessarily 



contains some unknown parameters. Next ip{q) should be plugged in the eq.( p.l| ) and the 



resulting expression should be minimized with respect to the parameters. This procedure 
gives the free energies of different mesophases as functions of temperature and composition. 



The molecular properties enter via the coefficients in eq. (|2.1| ). For the given temperature 
(below the ODT) and composition / the mesophase with the lowest free energy will be the 
stable one. The phase diagram is drawn by varying temperature and / over some ranges. 

The way to ascribe a given symmetry to the order parameter is to expand it in the full 
orthonormal set of functions with the given symmetry. The most popular choice is to use 
an expansion in plane waves: 

V'(q) = Y.A, (e^^-5q,k. + e-^^-5q,_k.) (3.3) 

i 

where {kj} is a set of vectors defining a lattice of the given symmetry, 0j and Ai are the 
arbitrary phases and amplitudes; 5ij denotes the Kronecker delta symbol. If the rotational 
symmetry in the system is not broken we can rewrite (|3.3|) in the following way: 



^(q) = E «^ E W'^'k k« + ^"'"'"^q -k«) (3-4) 

where the first sum runs over the Brillouin zones or coordination spheres in the reciprocal 
space; is the number of vectors in the i-th zone. In practice we consider only the first few 
terms in the series ( p.4|) . The mathematical reason for that is to simplify the minimization 
procedure. From the physical point of view, immediately below the transition temperature 
the appearing mesophase is not well developed and we do not need too many details in 
the representation ( p.4|) . (As a rule, the length of the g- vectors increases with the number 
of the zone. Therefore, the more terms are left in (|3.4]) , the more detailed density profile 
we get.) The next question is how many terms we should leave in ( |3.4| ). Unfortunately, 
there is no quantitative criterion where to cut the density expansion. In the past, several 
possibilities were considered. In the classical work |Q the expansion ( |3.4| ) was restricted 
to z = 1. The so-called second- harmonic approximation {i = 2) was formulated in 
Z^,[7| J25|J2E|JIP|] in order to study complicated bicontinuous phases. Up to several hundred 



harmonics were used in the numerical method of Matsen [§| (however, the authors studied 



the intermediate/strong segregation regimes and, therefore, our discussion is not directly 
apphcable to their model). Comparison of their results shows that higher harmonics not only 
result in some non-significant shift of the phase boundaries, but also stabilize or destabilize 
several morphologies. In such a situation the number of harmonics is crucial. 

In this paper we adopt the criterion derived in [II^,^. Instead of the mean-field theory 
Hamley and Podneks studied system with composition fluctuations. In their model the cut- 
off value for g„ (the g-vector from the nth shell) appears naturally as an adaptation of the 
Ginzburg criterion: 



-1 



< 1 (3.5) 



We apply this criterion to our mean-field theory in order to keep uniformity with the fluctu- 
ation approach (a similar criterion was used in [^). There are several ways to support our 
choice. One type of argument is connected with the applicability region for our theory. In 
the RPA there is essentially one lengthscale associated with the gyration radius of molecules. 
It is intrinsically connected with the typical wavevector q^... Therefore, the representation 
of the density profile by higher harmonics with g ^ g^, is in strong contradiction with the 



used theory. (A similar idea was formulated in a different way in p8[.) One can also specu- 
late that near the spinodal the inequality ( p.5| ) will always hold in the proper temperature 
range. The questions of applicability of ( p.5[ ) and its influence on the phase diagram will be 
addressed in future studies. 

The cut-off (|3.5| ) implies that only one harmonic should be used for LAM, HEX, BCC, 
SC and SL phases, while FCC, GYR and OBDD phases are represented by two-harmonic 
density profile. The basis vectors for LAM, HEX, BCC, SC and SL phases, as well as 
the first-shell vectors for FCC and OBDD (which coincides with the BCC first shell), are 
presented in PJl8[]. For the rest we use (more detailed information can be found in p9|): 
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FCC: fcf) = ^(1,0,0) 



OBDD: A;f) = ^(1,1,1) 



GYR: fcW = ^(l,_l,_2) 

fcf) = %f_l,_2,-l. 



V6 



4'^ = ^(1,1, -2) 
^lo"* — "^(-'^' ~2, — 1) 



1,(2) 


= ^(0,1,0) 


4'^ = ^(o,o,i) 


l(2) _ 




= ^(-1,1,-1) 


iL(2) 






_ 

At2 — 


^(-l'-l'-2) 


A:i'^ = ^(1,2,-1) 


_ 

^5 - 


= ^(-2>l>-l) 


^6 = :^(~2, —1, — 1 


^.(1) _ 
/tg - 




4'^ = ^(-1,2,-1) 


'^11 ■ 




A:« = ^(2,-l,-l) 



fcf ) = ^(0, 1, 1) fcf = ^(0, 1, -1) fcf = ^(1, 1, 0) 



A^f = ^(1, -1, 0) # = ^(1, 0, 1) A;f = 0, -1) 

Once we have defined the density profile (|3.4| ), we plug it into the equation for the free 
energy The integration is trivial since ipi^) is a sum of delta-functions. The resulting 

expression contains terms proportional to the vertex functions T3 and r4 which arguments 
are all possible combinations of two, three and four fcj*^ that add up to zero. For the phases 
represented by the two-harmonic density profile some of these terms will be the functions of 
the wave-vectors from the different shells. However, while calculating the vertex functions 
(Section 2 and Appendix A and B) we assumed everywhere that all wave-vectors have the 
same modulus g*. Correspondingly, we obtained the m-th order vertex function in the form 

rm(5'*ni, . . . , g^Kiim) = Tmix^:, h^"^^) (rij = — ; /i*^™^ and x are defined in the Appendix A). 

Qi 

This function cannot be directly used for the arguments with the different moduli. However, 
the following trick helps to avoid this problem. 

The criterion (|3.5| ) implies that qj = q^^ + Aj, where Aj ^ g*. Therefore, to the first order 
in A we can approximate the m-th order vertex function of q-s with the different length by 

m 
a=l 

where 

. . . , qa, . . . , g*n^ 



dqa 
10 



qa=q* 



Any vertex function is a symmetric function of its arguments. Therefore, inde- 
pendent on a. Let us now put all Aq, = A. This gives 

rm(qi, . . . , q^) ~ T^{q^wi, g^n^) + uikA (3.6) 

But this expression should be identical to: 

6x (3.7) 

where x = (g* + A)'^Rq = x=k + 6x, 6x ~ 2x* — . Comparing ( |3.71 ) with ( p.6| ) gives 

q* 

2x, dTmix^h^""^] 

K 



ax 



Finally, 

r^(qi, . . . , q„) ^ f ^(5: ^, h^-y) = r„(x., /^(-)) + ^^iziii^lZL^ 



(3.8) 



The angle variables h^^^ are calculated from the vectors g^Kiij. Here we assumed that the 
latter sum up to zero, which is only approximately true. The approximation (|3.8|) allows us 
to take into account the small differences in the length of the wave-vectors from the first 
and the second shells. 

The last step involves the minimization of the free energies with respect to the arbitrary 
phases (p^j^ . The algebra is labouring but straightforward. Denoting ai by a, and a2 by b we 
obtain 



Ffcc = ^a'T^ix,) + Sb'T^ i^jx^j + Ua'bT, |^-== - ij + E^'a^ + ^2^-^54 ^ 421-^0^5 
FoBDD = Ga'Ux.) + 4b'T, (^x.) + 8a'U0) + ^Ef'^a' + + 2Sf a^ft^ (3.9) 

Fgyr = 12a2r2(x,) + 6b^r2 Qx,) + Sa^TsiO) - 12a%r3 (^-^ - 1^ 

+86=^f3 {2V3 - 3) + SEf^a' + ^Efh' + 2eP a'b + 2EPa%' 

where the coefficients E are given in Appendix C. 
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The one-harmonic free energies for the LAM, HEX, BCC, SC and SL phases are presented 
in P,p!8|. As it was noticed in we have two choices for the FCC free energy, which 
arise from the different solutions of the phase minimization equations. The first solution is 
presented in (|3T9|). The second (NFCC in the notation of |^) was never found to be stable 
in our calculations. 

We want to emphasize that the second harmonic approximation was used extensively by 
many authors Ill8| , p4| - p^ , |l9| ,f^ . However, almost all of them used the two-harmonic repre- 
sentation of the density profile for all phases. On contrary, we apply the criterion (|0|) and 
represent only the FCC, OBDD and GYR morphologies in the two-harmonic approximation. 
This approach was also tried in several studies [p!0tp|,p7|. The main difference with our work 
is that we not only keep the complicated angle dependence of the vertex functions, but also 
take into account their dependence on the wave vectors moduli. 



IV. RESULTS AND DISCUSSION 

In this section we apply the developed formalism to particular systems. As an example 
we will treat dendrimeric tree-like structures defined as follows. The origin of the molecule 
is a branching point with ni A-blocks of the same length fi attached to it. All rii A- 
blocks are topologically equivalent. They form a so-called first generation of the molecule. 
The end of each block is also a branching point with n2 B-blocks of the same length /2 
attached to it. All ^1^2 B-blocks are topologically equivalent and form the second generation. 
The third generation again consists of A-blocks and so on till the g-th generation - the 
highest generation of the molecule. We introduce a (7- dimensional topology vector n = 
{ill, n2, . . . , Tig}. The i-th element of this vector defines the number of blocks (A-blocks if i 
is odd, and B otherwise) originated from eac/i block in the (z — l)-th generation. The lengths 
of the blocks within one generation are necessarily the same but can be different in different 
generations. It is therefore useful to measure the lengths of the A-blocks in terms of the 
length of the first generation = fi {fs = /2 for the B-blocks correspondingly). Thus, we 
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introduce the following (^-dimensional vector r = {1, 1, rs, . . . , r^}. The length of the i-th 
generation is given by 

fi = n{fAe^ + fB{l-e,)) (4.1) 

where 

{1 i is odd 
otherwise 

If we now introduce the total content of the A-component in the molecule p, then the 
following equations hold: 

9 i 
P = IaYI ^i^i n 

i=l j=l 

l-P = /BE(l-e.)r.nn,- (4.2) 
i=i j=i 

and for the given topology (fixed g, n and r) and the overall composition p, /a and fs are 
determined in a unique way. 

Next we want to be able to calculate the distance between two given blocks. We introduce 
the system of coordinates and number all blocks originated from some block on the (z — l)-th 
generation from 1 to n^. Then the position of some block is given by a set of numbers, that 
define the trajectory from the origin of the molecule to this block, and will be denoted by 
a (yf-dimensional vector R. For the i-th generation block all components of R higher than i 
are equal to zero, (see Fig.|l|). The shortest distance between the ends of two blocks with 
positions R*^^-* and R*^^^ is given by 

A = (/^e, + fsil - e.)) (l - ^H(^)-Hf),o) (2 " K^\o ' K<^\o) (4-3) 

We reformulate the framework of Sections 1 & 2 in the form of algorithm, depicted 
schematically as follows: 

1) Input: g, n and r 

2) Vary p and the ranges of interest and for each point: 
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- calculate the vertex functions 

- calculate the free energies for the mesophases 

- find the stable mesophase (the one with the lowest free energy) 

3) Output: the phase diagram 

The calculation of the vertex functions is preformed in the fashion described in Section 
1. For example, to calculate the third order correlation function Gaab we need to move two 
labels over all A-blocks, and one label over all B-blocks. For each configuration of the labels 
we need to determine the type of diagram and then use the length of the blocks with the 
labels for f^, //3, • • • and the distances between those blocks for Ai, A2, ... in the equations 
from Appendix B. The resulting number should be added to Gaab- 

To distinguish among the different types of diagrams (line, fork, star) we introduce the 
concept of pathways. A pathway from one block to another one consists of the positions of 
all blocks lying in-between those two (so, the pathway is a set of maximally 2{g — l) vectors). 
Now suppose we know that the labels (three or four of them) occupy three different blocks. 
According to Appendix A they can form a line- or a star-type diagram. In the first case one 
of the blocks belongs to the pathway between the others, while in the second case it does 
not. Checking this condition allows us to distinguish among different types of diagrams (In 
the 4-block case the number of conditions increases but the general idea remains the same). 

The procedure of the free energy minimization with respect to the amplitudes a and b is 
the same as that of Leibler and is not presented here. We proceed with the results for 
different dendrimers. 

Figs.0-^ show the phase behaviour of various architectures. Despite the differences in 
topology, number of blocks and branching points, all systems behave in a similar way. 



Their phase diagrams are very close to those of the diblock M, triblock [11], stars and 



graft copolymers ||T2|,|T3[. The total number of blocks and branching points influences the 



spinodal temperature, while the ratio ua/ub and the number of branching points influences 



the asymmetry of phase diagrams. (A similar observation was made in |jT2| where the authors 
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connected the increase of complexity of molecules with the shift of the spinodal temperature.) 
Moreover, the relative sizes of the stability regions for different phases are almost the same 
for all systems, while the types of phases and their order is completely independent of 
the topology of molecules. Therefore, we conclude that all mean-field model systems with 
Gaussian chains will have phase diagrams similar to Figs]^-|^. This is in agreement with 
the recently proposed constituting block copolymer hypothesis [^] which explains the phase 
behaviour of complicated molecules on the basis of the behaviour of their constitutive blocks. 
For example, in the case of the two-generation dendrimer from Fig.^], the A1B2 miktoarm- 
stars are the building blocks of the molecules and their behaviour qualitatively corresponds 
to the phase diagram of the {2, 2}-dendrimer. Since we can always construct a complicated 
molecule from simple ones, like diblock, triblock and so on, the resulting phase behaviour 
is going to be typical. However, unlike the general theory of Sections 2 and 3, our model 
system does not consider structures where the ends of the blocks are the branching points 
with both A and B blocks attached to them. (This is necessary, for example, to model comb 
copolymers.) Therefore we cannot be sure that such systems will also have typical phase 
diagrams. There is also evidence of the opposite: Dobrynin and Erukhimovich found that 
the theory fails for the homo-stars with n > 5 and the corresponding phase diagrams 

contain regions where the phase behaviour cannot be predicted within the used approach 
|r3| . Nevertheless, the other phase diagrams show the same typical behaviour. 

We also want to notice the complete absence of the OBDD morphology. This is not 
surprising even for the {4, 4, 4}-dendrimers with the 4-fold symmetry of molecules. As it 
was shown by Matsen and Bates 0, the OBDD phase has larger interface in comparison 
with the GYR phase and could be stable only at very low temperatures. It is however 
screened by other morphologies (BCC and HEX). 

Another interesting feature is a notably big region of stability of the bicontinuous gyroid 
phase in Fig.^. For such molecules the highly branched structure and the relative confor- 
mational freedom of the end A-B blocks favor the GYR phase to the LAM phase. This is 
not a feature of the {3, 3, 3}- and {4, 4, 4}-dendrimers, for example, because there the outer 
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shell of the molecule consists only of the A blocks. As the AB dendrimers were not studied 
experimentally yet, we recommend the systems like {5, 1, 3, 1} to produce relatively easy the 
bicontinuous gyroid phase. 

At the end we want to mention that our tree-like model is capable of describing stars, 
polyblocks, umbrellas, and many other architectures. Because of space limitations those 
results are not presented here. 



V. CONCLUSION 

In the present paper we have shown how to calculate phase diagrams for systems with 
arbitrary topology. The developed method is a generalization of the work by Leibler [Q. 
The monomeric correlation functions are calculated with the help of a special diagrammatic 
technique which allows us to treat an arbitrary molecular topology. Apart from the assump- 
tion that the polymers obey Gaussian statistics even below the order- disorder transition 
temperature, the calculation of the vertex functions is approximation-free. The free ener- 
gies for different mesophases are calculated within the second-harmonic approximation as it 
was formulated by Hamley and Podneks . The influence of this approximation and its 



relationship with the random-phase approximation requires additional study. 

The architecture of polymer molecules is known to influence the phase diagram only 
slightly This was once again confirmed by our studies of highly-branched model 

dendrimers. We found that the topology of molecules influences the spinodal tempera- 
ture and asymmetry of the phase diagrams, but not the types of phases and their or- 
der. The typical phase diagrams are presented in the Results and Discussion. That 
model system, however, does not cover all possible architectures. Many more phase di- 
agrams together with the C-code used to obtain them can be found on the Internet: 
http://rugmd4 ■ chem. rug, nl/ ~morozov/ research, htrrn . 



Our method can be a starting point for a fluctuation theory of systems with arbitrary 
topology. This will require some straightforward generalization of the Hartree approximation 
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by Brazovskii p2| , |33| and Fredrickson and Helfand p4|, since for some systems (see, for 



example, Figs. 4 and 5) the 4-th order vertex function r4^{hi,h2) exhibits a strong angle- 
dependence and, therefore, cannot be approximated by the constant value r4(0,0) or any 
other. One of the possible ways for such a generalization was proposed by Dobrynin and 
Erukhimovich |35[ |. 

VI. ACKNOWLEDGMENTS 

We thank A. Zvelindovsky, M. Reenders and I. Reviakine for reading the manuscript. 



17 



APPENDIX A: DIAGRAMMATIC TECHNIQUE 



We consider a polymer molecule consisting of ua blocks of type A and hb blocks of 
type B. The total number of monomers in the molecule is A^. We shall use Greek letters to 
number the blocks: a, P,j,6 = 1 . . .ua + ub- The length of the a-th block is La = faN. 
We use solid lines to denote blocks and dashed lines to denote arbitrary distance between 
blocks. The circles denote the labeled monomers. We use s, s', s" . . . for the positions of 
the labels on their blocks. A diagram with n labels will be called the n-th order diagram. 
The labels act as sources of momentum. A circle labeled with a,b,c. . . emits q^, q^, qc . . . 
respectively. The total momentum of any diagram should be zero, therefore all external 
lines should carry zero momentum. We also introduce a convention (which does not affect 
our results) that all momenta flow from left to right along the chain. This can be illustrated 
by the following: 

qi qi + q2 o 



qi 



q2 



qs 



A given number of blocks can form topologically different diagrams: 



2_ line ^ line star j_ 



J _ line fork \ \ star \ 

. _ _ _ _ _ 



Necessarily, all momenta flowing in and out of each branching point should add up to 
zero. 

Following 0] we assume that all wavevectors have the same length |q| = and introduce 
the following angle variables: 
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for the 3rd order diagrams: = <llhf^ , h^^^ — (l,l,h) 
for the 4th order diagrams: (q^ + q^)^ = qlh[f 



where 



4 


hi 


4 - /li - /l2 


/l2 




4 




4-/ii 


hi — /i2 


h2 


4 




/i2 


4 - /ii - /l2 


hi 


4 



The matrix h^^^ has two important symmetry properties: h\f — h^^^ and h^^^ — ^^^^ 



for 



'i ^ 3 1^ k ^ I. These properties will be extensively used in Appendix B. 
For any internal line we subscribe the Gaussian propagator: 



V{1) 



-—I 

e 



where Rq is the radius of gyration of an ideal chain of N monomers. The product of the 
propagators is integrated over all possible positions of all labels on their blocks. The result 
should be then multiplied with the symmetry factor arising from all possible permutations 
of labels that result in topologically the same diagram and keep labels on the blocks they 
initially occupied. Finally, the normalization constant 1/N is added. 

The following functions will appear as the building blocks for all diagrams: 

1 - e-f'' 



Ff 
Tfih) 

J Ah) 
Vf{h) 
Ti{f,h) 



X 



fh 



h 



X^ 



X^ 



1 



7 

-fx 



f 



-hfx 



h h-1 h-h'^ 

hfx + e-^f"" - 1 , 



1 - e-f^ - f. 



xe 



-fx 



h^ 



fx 



3- fx 



fx' 



h ^ h 



: fx 



+ e 



-fx 



2h-3 



-fhx 



fx + e"^^ - 1 



2h + l 
~h?~ 



fx' _ 



fx' 



e-^^ + + /x - 3 + 3e--^^ 



+ 3 - - (3 + /x)e--^^ 
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T2{f,h) 



1 



-fhx 



h-2 



h h{h-iy {h-iy 

fx' 



fx 



1-e 



-fx 



fxe 



-fx 



h-1 

I 2e-/^ + /xe"^^ + fx-2 



All functions of the angle variable h are defined in the following way: 

...(h) = [ h^O,l I h = l I h = ] 

Some of these functions were previously defined in P. |TT|j36[| . Separate elements of the pre- 
sented diagrammatic technique were used in [p^ , p^^ -|39| • 



APPENDIX B: CALCULATION OF DIAGRAMS 



+ 



s' 7VA s 



— 



2 



s" s' NA s 



-•— + • • + 7— • -? 



s" NAi s' NA2 s 



'AT A3 
I 

s NAi , NA2 

+ Tr-« ^ - -, 



II) 



La 



^ E J/d^l ds'l ds"V{hf{s - s') + hf\s' - s")) 
= 2N' [2g2{U,h)+g2{fa,l)] 



— J ds I ds' I ds" 



JO 



V (/if {s' - s") + [Lp -s' + Nd. + s] 
+ V {hf\s' - s") + /i(3)(L^ _ + iVA + s)) 
N'e-^'''-'^:FfShf)\jf,{h) + Jf^{hf) 
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^ ds r ds' r ds" V [hf (s" + iVAi + s') + hf {Lp -s' + NA2 + s) 
i\ Jo Jo Jo ^ 



^ ds j^' ds' ds" V [hf {s + iVAi) + /if {s' + 7VA2) + hf {s" + TV A3)) 



+ 7 



+ 



b*->a 



s'" NAi s" NA2 

+ ^ ;r-» 



,NA3 



s'" NAi s" s' NA2 



+ TT-* - „ > • • + • _ ± - 

^ b "a "a ATAi NA2 1^ b 



s'" NAi s" NA2 s' NAa s 
+ -T-* T-* 3-* TT*- 



NA4, 



s'" NAi s" NA2 



r-* 



+ 



'AT A3 



WAi 



ArA2 ' 



AfAs 



iVA4 



'iVAs 



r: 




r(s"-s"' + hf^{s'-s")+s-s') 
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4 E / ds ds' ds" / ds'" 
N ^J^^^Jo Jo Jo Jo 

p (s" _ s'" + h'^Jis' - s") + L^-s' + NA + s) 



4 rl^a fS l-La fs" 

i- ds ds' ds" / ds'" 
N Jo Jo Jo Jo 

p _ s'" + {Lp - s" + 7VA + s') + s - s') 

I- 1"^" ds r ds' f^^'ds" ['"'ds'" 

N Jo Jo Jo Jo 

V (s'" + TVAi + s" + /iiJ (L^ - s" + iVA2 + s') + s - s') 



4 E / ds ds' ds" / ds'" 

V [s'" + iVAi + s" + - s") + L)3-s' + NA2 + 



/•ivcK f'^f3 r^'v c^" 
4t ds ds' ds" / ds" 
N Jo Jo Jo Jo 

V(s + iVAi + s' + NA2 + h^^i (iVAs + s'") + s" - s'") 

1 rLoL rLg i-L^ rLs 

— ds ds' ds" / ds'" 
N Jo Jo Jo Jo 

V [s'" + A^Ai + s" + /iit^(L^ - s" + A^A2 + s') + - s' + NA^ + 

1 fLa fLfj rLj rLg 

i- ds ds' ds" / ds'" 
N Jo Jo Jo Jo 

V{s + NA^ + s' + iVAs + s'" + NAi + s" + /liJ (L^ - s" + NA2 

^3^-.(A,+.^)A.+A3+A,)^^^^^^^/,.^^^(^(4) _ 1)^^^ 
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1 rLa fLg fL^ fLg 

IX) ^4- ds ds' / ds" / ds'" 
N Jo Jo Jo Jo 

V[s + NA^ + s' + NAs + /iijATAs + s" + NA2 + s'" + A^Ai) 

= F, F, F, Ff p-^(Ai+A2+A3+A4+/iWA5) 



APPENDIX C: COEFFICIENTS S 



~ ( Aa \ ~ ( ^ Aa \ 

We use the following notation: ^ — ^1, /i2 = ^4 ^ — h^'^^ . 

\a=l ^* / \a=l ^* / 

FCC: 

SS^) = 6f 4 (0, 0, ^) + f4 (0, 0, 0) + 2f 4 (0, ^, ^) 
Hf^=4f4 (4g,0, 2) + f4 (4g,o, 0) 

) = 3f 4 (2Q, 0, ^) + f 4 (2g, 2, ^) + f 4 (2g, ;2) + f 4 (2Q, I ;2 

OBDD: 

) = 2f4 (0, 0, 2) + 8f 4 (0, 0, 1) + f 4 (0, 0, 0) + 414 (0, 1, 1) 
) = f 4 (4Q, 0, 0) + 2f 4 (4Q, ^, ^) + 6f4 (4g, 0, ^) 

) = 4f4 (2g, 0, 2) + 6f4 (2g, 0, z) - 2f4 (2^, 2, ^) - 6f4 (2g, 2, ^) + f4 (2g, 

+f 4 (2g, 4, ^) + f 4 (2g, 2f 4 (2g, ^, ^ 



GYR: 



Hp^ = f 4 (0, 0, 0) + 4f 4 (0, 0, + 4f 4 (0, 0, 1) + Sf 4 (o, 0, ^) - 4f 4 (o, i 
- (°' % D + (°' °' t) + (°' I D + (°' °' D 

sf ' = 8f4 (4g, 0, 1) + 2f4 (4g, 0, 2) + 4f4 (4g, 1, 1) + f4 (4g, O, 0) 



1 2 
3 
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f ' = 2f4 (q,0, - 2h (q,0, 5) + fi (q, i, i) + tt (q, |, I) + 2f4 (o, 1, 
- 6f 4 (q, 2, i) - 2f 4 (q. 2. ij") - 2f 4 (Q, 2, v^) - 2f 4 (q. i. i) 



4 



j +2f4 (^Q,V3,2--^j 

- 2f4 (g, I ^ + v^) - 2f 4 (g, + 2f 4 (g, y , 2 - x/a) 

-2f4 (^Q,2-VS,-^^ +2f4 (^g,-^,2-^j +f4 (g,2V3,2- V3) 

= 6f 4 (2g, 0, 2) - 2f 4 (2Q, 2, -^^ - 2f4 (2g, 2, v^) - 2f4 (2Q, ^, 2 - i= 

- 2f4 (2g, I 2 - ^) - 2f4 (^2g, 1, 2 + j - 2f4 (^2Q, ^> 2 + -^j 
+ 12f4 |^2g, 0, 2 + j + 6f4 |^2g, 0, 2 + - 2f4 (2g, 1,1 + 

- 2f 4 (2g, 1 ^ + x/s) + 2f 4 (2g, 3, \/3 - 1) + 2f 4 (2g, 4, 2 - \/3) 

- 4f 4 |^2g, 2 - \/3, j - 4f4 |^2g, 2-V3, j + 2f4 (2g, 2-V3, 2V3) 

+ 12f 4 (2g, 0, 2 - v^) 
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FIGURES 




a) b) 

FIG. 1. Model dendrimeric tree-like structure. The empty circle denotes the origin of the 
molecule. The filled circles denote the labeled blocks, a) Here we illustrate the introduced system 
of coordinates. The positions of the labels A, B and C arc R(^) = {1,2,2}, R^-^) = {2,0,0} and 
R(*-^) = {2,3,2}, respectively, b) The bold lines denote the pathways between the labeled blocks. 
The pathways are: PW(A,B)={{1,2,0},{1,3,0}} and PW(C,D)={{2,2,0},{2,0,0},{3,0,0},{3,2,0}}. 
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FIG. 4. Phase diagram for the 5 = 3, n = {3, 3, 3} and r = {1, 1, 1} dendrimer. 
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FIG. 5. Phase diagram for the = 3, n = {4,4,4} and r = {1, 1, 1} dendrimer. The cartoon 
of the dendrimer is not in scale: all lengths of the A-blocks should be the same. 
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